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ABSTRACT 

We describe the discovery of a likely brown dwarf (BD) companion with a minimum mass of 31.7 
± 2.0 Mj up to CSC 03546-01452 from the MARVELS radial velocity survey, which we designate 
as MARVELS-6b. For reasonable priors, our analysis gives a probability of 72% that MARVELS- 
6b has a mass below the hydrogen-burning limit of 0.072 M©, and thus it is a high-confidence BD 
companion. It has a moderately long orbital period of 47.8929^q'qq62 days with a low eccentricty of 
0.1442lg;Qg 7 3, and a semi-amplitude of 1644^13 m s _1 . Moderate resolution spectroscopy of the host 
star has determined the following parameters: T e g = 5598 ± 63, logg = 4.44 ± 0.17, and [Fe/H] = 
+0.40 ± 0.09. Based upon these measurements, GSC 03546-01452 has a probable mass and radius of 
ill* = 1.11 ± 0.11 Mq and R* = 1.06 ± 0.23 Rq with an age consistent with less than ^6 Gyr at a 
distance of 219 ± 21 pc from the Sun. Although MARVELS-6b is not observed to transit, we cannot 
definitively rule out a transiting configuration based on our observations. There is a visual companion 
detected with Lucky Imaging at 7.7" from the host star, but our analysis shows that it is not bound 
to this system. The minimum mass of MARVELS-6b exists at the minimum of the mass functions for 
both stars and planets, making this a rare object even compared to other BDs. It also exists in an 
undcrdense region in both period/eccentricity and metallicity /eccentricity space. 
Keywords: stars: individual (GSC 03546-01452) 
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1. INTRODUCTION 

Radial velocity (RV) surveys have provided a wealth of 
exoplanet discoveries aroun d sun-like stars i n rece nt years 
(California Planet Survey, iHoward et al.l (120101 ) Lick - 
Carnegie Exoplane t Survey, Hagh ighipo ur et al.l (|2010D : 
CORA LIE survey. lUdrv et all (120001): and the HARPS 
survey, iMavor et all ( 20031 ) to name a few), but they 
have not found a correspo ndingly larg e number of brown 
dwarf (BD) companions (jReid fe Metchevl I2008D . BD 



companions lie on the mass spectrum between planets 
and stars and are defined as being between l3Mj up and 
75.5 M,) „„ (based on the deuterium and hydrogen fusion 
limits) (jChabrier et al.l[2000t ISpiegel et al.ll20lTI ). The 
lack of BD s within 3 AU of their host star was first rec- 
ognized in iMarcv fc Butlerl (|2000( ). and is known as the 
BD "desert". This result is unlikely to be due to obser- 
vational bias because the RV semiamplitudes of BD arc 
many hundreds to a few thousand meters per second, 
which is easily detectable by these surveys (|Patel et all 
120071 ). Thus the lack of BD at close to moderate dis- 
tances from their respective host stars points to an ex- 
planation based in BD formation mechanisms. 

Although the formation of low- mass companions (plan- 
ets through low mass stars) to sun-like stars is still an 
area of active research, in overview there are two main 
mechanisms: planets form from a protoplanetary disks 
and stellar companions [in a similar range of separations] 
from molecular cloud fragmentation. Given the mass 
range of BD companions they could form through cither 
mechanism (or both). Understanding the origin of the 
BD desert can put major constraints on the upper mass 
limit for companion formation in protoplanetary disks, 
and a lower mass limit on formation via fragmentation. 

Recent efforts to quantify the frequency of companions 
as a function of mass in the BD desert have found that 
the overall frequency of BD companions at close to mod- 
erate d istances from their host star ( < 10 AU) is less than 
< 1% dGrether fc Lineweaverl 120061 ). and more recently 
0.6% (jSahlmann et all 120111 ). This val ue is low com- 
pared to ~ 7% for planetary companions (jUdrv fc Santos! 
120071) and ~ 13% for ste llar companions in a sim- 
ilar range of separati on s dDuauennov fc Mavorl 119911 : 
iHalbwachset all l2003f ). iGrether fc Lineweaverl (J2006D 
went a step further and defined the driest part of the 
desert to be where there was a minimum in the number 
of companions per unit interval in log mass; they found 
this position to be at a companion mass of 31^-jjMjnp. 

This is the fifth paper in this scries looking at low-mass 
companions to sun-like stars from the th ird generation of 
the Sl oan Digital Sky Survey (SDSS-III; lEisenstein et all 
(J2011I )) Multi-Object APO Radi al Velocity Exoplanet 
Large- Area Survey (M ARVELS; IGe et al.lBool. I200J 
iGe fc Eisenstem|[2009l ). The primary goal of this series 
is to provide a detailed set of well-characterized compan- 
ions with minimum masses and separations in or near the 
BD desert, which can be used by future meta-analyzes. 
Ultimately, this is the sort of groundwork that must be 
done in order to understand the extent and aridity of the 
BD desert. 

The MARVELS survey measured radial velocities of 
3,300 unique FGK type stars. MARVELS is a large 
survey looking for RV companions around bright stars 
(7.6 < V < 12) with periods below 2 years with well 
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Figure 1. Top: The original detection Lomb-Scargle pcriodogram 
for GSC 03546-01452 based on the MARVELS data. The two RV 
points with the lowest spectrum flux were removed for this initial 
pcriodogram (as marked in Table [TJ. The peak in the pcriodogram 
(marked by the red dashed line) is quite significant and is at 47.842 
days, close to the final adopted value of 47.893 days for the com- 
bined MARVELS and TNG/SARG data. Bottom: A combined 
pcriodogram for all 60 stars on the MARVELS plate. The highest 
power point in each frequency bin was removed, and the remaining 
were averaged. The red dashed line shows the peak period from 
the periodogram on the left. There is no obvious systematic peak 
on the plate periodogram at the peak period of the GSC 03546- 
01452 periodogram. 

characterized biases; see ILee et all (1201 ID for a descrip- 
tion of the survey desig n. O ther papers in th is se- 
ries (iFleming et all 120101 I20H iWisniewski et~aTl l20ll 
iMa et al.ll2013al Mack et al. 2013, in press; Jiang et al. 
2013, submitted) have helped fill in our understanding 
of the BD desert, and provided warnings to some of the 
pitfalls inherent in these analyzes. 

We will discuss observations of the star GSC 03546- 
01452, which has a companion with a period of ~ 47 days 
and with a minimum mass of 31.7 ± 2.0Mj up placing it 
near the most "arid" region of the BD desert. In Section 
2 we discuss the photometric and spectroscopic observa- 
tions and basic data processing. In Section 3 we discuss 
the analysis of this data. Section 4 contains a discussion 
of the results and places MARVELS-6b within the larger 
context. Section 5 summarizes our results. 

2. OBSERVATIONS AND DATA REDUCTION 



MARVELS-6B 



The initial detection and orbitial characterization of 
MARVELS-6b came from the MARVELS Survey RV 
data. Once it became a strong candidate further RV, 
photometric, time-series photometric, and spectroscopic 
data were taken to confirm this candidate. Each of these 
data sets will be discussed below. 

2.1. Initial Identification of MARVELS Candidates 
2.1.1. Survey Summary 

MARVELS is a multi-epoch radial velocity survey de- 
signed to detect radial velocity companions around FGK 
type stars in a magnitude range of (7.6 < V < 12). 
It uses a disp ersed fixed-delay interferom eter (DFDI; 
Ge et al.ll2009h on the SDSS 2.5m telescope (jGunn et all 



20061) . The DFDI method was i ntrod uced for use in 
a multi-object RV survey by [Gl (|2002l ). A single ob- 



ject version of a DFDI instrument was success fully used 
to de tect a hot Jupiter around HD 102195 (jGe et al.l 
120061). The DFDI instrument principle was described by 
(lGell2002HGe et al.l[200l IErskindl200l Ivan Evken et all 
120101: iWanget allboill ). The MARVELS interferom- 
eter de lay calibrations were described in iWang et all 
(l2012aibH 

Each MARVELS observation consists of 60 stars 
spread across 120 spectra (2 spectra for each star). The 
MARVELS survey started in Oct 2008 and ran through 
July 2012. It was divided up into two observing sets: 
year 1-2 fields which were observed from Oct 2008 till 
Jan 2011 and year 3-4 fields which were observed from 
Jan 2011 till July 2012. 

The year 1-2 data set con sisted of 43 unique fields (11 
of which are in the Kepler [Borucki et al.ll20T(il |) field of 
view. Year 3-4 field set contained 13 fields (with 1 field 
overlapping the year 1-2 field set.) This leads to a total 
of 3,300 unique stars with > 18 epochs. 

2.1.2. Radial Velocity Analysis software 

The large number of RV targets required development 
a software package to easily display and characterize the 
radial velocity curves from the MARVELS survey. One 
of the primary tools used in this process is an IDL-based 
Keplerian model fitter known as MPRVFIT3 MPRV- 
FIT starts with the input of RV data in the form of Ju- 
lian date, RV, and RV error (the Julian date and RVs 
can be of any type that is appropriate to the analy- 
sis). In the case of MARVELS, the dat a from the two 
beams have been combined as described in lFleming et all 
(|2010l) . The software then applie s a modified version 
of th e Lomb-Scar gle periodogram (jLombl 119761 ; IScargld 
I1982D described in ICummingl ( 20041 ). False alarm proba- 
bilities wer e assigned to t he highest peaks using the for- 
mulation of lBaluevI (|2008l ). Once a number of high prob- 
ability peaks are identified, the frequency space between 
the peak and the next nearest frequency point in the 
periodogram is sub-divided into 10 frequency steps (on 
both sides of the peak) . A Kepl erian model is then fit to 
those frequencies using MPFIT (|Markwardtll2009l) . MP- 
FIT is a Levenberg-Marquardt non-linear least squares 
fitter implemented in IDL. The x 2 statistic is determined 
for each fit, and the best chi-squared fit of the grid is re- 
tained. Finally, the best fit models are plotted with the 

24 http://www.vandcrbilt.edu/AnS/physics/vida/mprvfit.htm 





Table 1 




Summary of Radial Velocity 


Data 


BJD 


RV 


RV err 


Source 


(days) 


(m s- 1 ) 


(ms- 1 ) 




2454956.915966 


2039.20 


72.10 


MARVELS a 


2454957.946510 


2150.89 


43.11 


MARVELS 


2454958.912457 


2206.71 


39.27 


MARVELS 


2454959.927110 


2291.25 


36.06 


MARVELS 


2454962.924107 


2568.84 


36.54 


MARVELS 


2454963.890111 


2559.01 


48.00 


MARVELS 


2454964.918085 


2529.69 


39.97 


MARVELS 


2454965.910038 


2579.81 


43.36 


MARVELS 


2454984.835892 


-797.93 


41.72 


MARVELS 


2454986.892198 


-736.07 


43.95 


MARVELS 


2454988.905306 


-444.95 


40.44 


MARVELS 


2454990.799823 


-154.83 


64.34 


MARVELS 


2454993.804470 


255.99 


43.11 


MARVELS 


2454994.793741 


543.09 


49.32 


MARVELS 


2454995.898027 


594.44 


43.11 


MARVELS 


2455014.802962 


2438.46 


51.41 


MARVELS 


2455020.845356 


1639.58 


39.45 


MARVELS 


2455021.851131 


1381.28 


50.68 


MARVELS 


2455023.835931 


792.32 


53.52 


MARVELS 


2455024.788660 


463.02 


51.78 


MARVELS 


2455484.614289 


2165.40 


59.70 


MARVELS 


2455436.578251 


1901.69 


24.22 


TNG/SARG 


2455460.467627 


-841.33 


4.49 


TNG/SARG 


2455460.489293 


-833.73 


4.44 


TNG/SARG 


2455460.511619 


-836.81 


4.76 


TNG/SARG 


2455521.551045 


722.14 


88.25 


MARVELS b 


2455725.687910 


1962.38 


6.86 


TNG/SARG 


2455760.505697 


66.39 


9.99 


TNG/SARG 


2455760.694598 


98.71 


5.87 


TNG/SARG 


2455791.532669 


122.53 


5.96 


TNG/SARG 


2455791.583907 


82.56 


5.52 


TNG/SARG 


2455844.434368 


-834.69 


5.08 


TNG/SARG 



a Second lowest flux spectrum 
b The lowest flux spectrum 

data points for easy reference. For the MARVELS sur- 
vey, the two best fits for each target were displayed, and 
the candidates were chosen based on the folded and un- 
folded radial velocity curves, as well as the significance 
of the periodogram. An example of this periodogram for 
GSC 03546-01452 can be found in Figure [TJ 

2.2. Radial Velocities 

2.2.1. SDSS-II MARVELS Radial Velocities 

Differential RV observations of GSC 03546-01452 were 
acquired during the first two years of the SDSS-III MAR- 
VELS survey. A total of twenty-two observations were 
obtained over the course of 565 days. As discussed in 
Section l2~j~Tl the MARVELS survey uses the DFDI tech- 
nique which introduces an interferometer into the light 
path. As a result, each 50 minute observation includes 
two fringed spectra, or beams, one from each arm of the 
interferometer. Each spectrum spans a wavelength range 
roughly 500-570 nm with a R~12,000 resolution. Each 
beam is processed individually through the MARVELS 
pipelin e following the methods described in iLee et all 

poll . 

The formal errors derived from the MARVELS pipeline 
are known to be underestimates of the true error. This 
systematic underestimate can be partially corrected for 
by using the fact that 60 stars (120 beams) are taken 
in each observation. Following the method outlined in 
iFleming et all (|2010f ). it is expected that most stars in 
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an observation plate arc radial velocity stable, so the 
median RMS is a reasonable estimate of the systematic 
errors. A quality factor (QF) is derived for each beam, 
which is the ratio of the radial velocity RMS to the me- 
dian formal error bar, and the median QF is found for the 
plate. The formal errors arc multiplied by this QF (2.334 
for GSC 03546-01452) resulting in the error bars used for 
the analysis of these observations. During the RV analy- 
sis, discussed in Section [3~2l it was determined that these 
scaled error bars were themselves overestimated, leading 
to a reduction of factor of 0.5794. The final differential 
RV measurements and final scaled error bars (approxi- 
mately 40% larger than the formal uncertainties) for GSC 
03546-01452 are presented in Table [U 

2.2.2. TNG Differential Radial Velocities 

Once GSC 03546-01452 was determined to be a can- 
didate for additional investigation, spectroscopic ob- 
servations were conducted wi th the Spettrografo Alta 
Risoluzione Galileo (SARG) (jGratton et al. I I2001D on 
the 3.6m Tclcscopio Nazionale Galileo (TNG) located 
at Roque de Los Muchachos Observatory (ORM). All 
the spectra were acquired using the same instrumen- 
tal configuration: a slit with a sky-projected width of 
0.8" achieving a resolving power of R = 57000; a yellow 
cross-dispersing grism providing the wavelength range 
462 < A < 792 nm. 

A total of 15 spectra were taken (11 with and 4 without 
the iodine cell inserted in the light path). One spectrum 
with the iodine cell proved to have too low S/N, leaving 
only 10 spectra suitable for differential velocity measure- 
ments. The spectra were processed using the standard 
IRAFp^l Echelle reduction packages. The S/N per resolu- 
tion element at 550 nm ranges from 50 to 130. The spec- 
tra without 12 lines served as a reference for measuring 
the relative radial velocities of the 10 spectra with super- 
imposed 12 absorption lines. The technique adopted t o 
derive RV values is described in lMarcv fc Butlerl (J1992D . 
For details on ho w the technique was implemented on 
SARG spectra, see lFleming et al.l (|2012[ ). As is the case 
with the MARVELS RV data, the formal uncertainties 
for the TNG/SARG data were overestimated, and so the 
error bars were thus reduced by a factor of 0.2224. The 
final differential RV and reduced error bars are listed in 
Table [U 

2.2.3. TNG Absolute Radial Velocities 

All 15 spectra taken with the TNG/SARG instrument 
were used to calculate absolute RVs. The stellar spec- 
tra were cross-correlated with a high resolution solar 
spectruirPI. The cross-correlation was done using the 
SARG red CCD spectrum with a wavelength coverage of 
6200-8000 A. Due to the non-simultaneity of the ThAr 
wavelength calibration exposures it is possible the slit il- 
lumination varied between the science spectra and the 
calibration spectra, which could result in the RV mea- 
surement being affected by systematic errors. To par- 
tially mitigate this error source we calc ulated the cross- 
correlation function of the telluric lines (jGriffin fc Griffinl 

25 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

26 http://bass2000.obspm.fr/solar_spcct.php 



Table 2 

SARG/TNG Absolute Radial Velocities 



BJD 


RV 


RV Err 


days 


(km s- 1 ) 


(km s- 1 ) 


2455436.550647 


-12.63 


0.13 


2455436.578251 


-12.73 


0.13 


2455460.467627 


-15.43 


0.13 


2455460.489293 


-15.35 


0.13 


2455460.511619 


-15.41 


0.13 


2455698.630462 


-15.28 


0.13 


2455698.654826 


-15.30 


0.13 


2455725.687910 


-12.48 


0.13 


2455760.505697 


-14.39 


0.13 


2455760.694598 


-14.32 


0.13 


2455791.532669 


-14.36 


0.13 


2455791.583907 


-14.22 


0.13 


2455791.610041 


-14.43 


0.13 


2455791.635364 


-14.52 


0.13 


2455844.434368 


-15.48 


0.13 



119731 ) around 6900 A with a numerical mask. This exer- 
cise results in a RV correction of a few hundred m s _1 . 
This correction plus the barycentric correction were ap- 
plied to the radial velocities and the resulting values are 
shown in Table [2] 

To validate this method, observations of the RV stan- 
dard star HD3765 were used. The standard was observed 
at 8 epochs with a mean result of -63,155 m s _1 ± 128 
m s . This result is within a 1-er agreement with the 
mean 34 observations taken from the ELODIE archive of 
-63,286 ± 49 m s _1 . For these observations, the system- 
atica caused by the non-simultaneity of the ThAr wave- 
length calibration is the dominant error source, so we 
adopt ± 128 m s _1 for the error in these measurements. 

2.3. High- Resolution Spectrum for Stellar Classification 

Two R -31,500 optical (~3600-10,000 A) spectra of 
GSC 03546-01452 were obtained on UT 2010 June 20 
with the Apache Point Observatory 3 .5 m telescope and 
ARC Echelle Spectrograph ( ARCES: IWang et al.l 120031) 
to enable accurate characterization of stellar fundamen- 
tal parameters (T e s, log g, and [Fe/H]). The two spec- 
tra were obtained using the default 1.6 x 3.2 arcsec slit 
and an exposure time of 1200 s for each spectra (for 
a combined exposure time of 2400 s). A ThAr lamp 
exposure was obtained after each of these integrations 
to facilitate accurate wavelength calibration. The data 
were processed using standard IRAF techniques. Fol- 
lowing heliocentric velocity corrections, each order was 
continuum-normalized, and the resultant continuum nor- 
malized data from each observation were averaged. The 
final spectrum yielded a signal-to-noise ratio (S/N) of 
approximately 110 per resolution element in the region 
around 6000 A. 

2.4. Photometry 
2.4.1. HAO Absolute Photometry 

We used the Hereford Arizona Observatory (HAO), 
a private facility in southern Arizona (observatory code 
G95 in the IAU Minor Planet Center) , to measure multi- 
band, absolute photometry of GSC 03546-01452. HAO 
employs a 14-inch Meade Schmidt-Casscgrain (model 
LX200GPS) telescope, fork-mounted on an equatorial 
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wedge located in a dome. The telescope's CCD is an 
SBIG ST-10XME with a KAF-3200ME detector. A 
10-position filter wheel accommodates SDSS and John- 
son/Cousins filter sets. GSC 03546-01452 and standard 
stars were o bserved with Johnson B and SDSS u', g', 
r', i' filters (jFukugita et al.l 119961 ). For the Johnson- 
Kron-Cousins b ands, standard stars are ta ken from the 
list pu blished bv lLandolt fc Uomotol (|2007l) and lLandoltl 
(|2009() . For the SDSS ba nds, standard stars are taken 
from the list published by iSmith et al.l <|2002f ) . Observa- 
tions were conducted on three dates: 24 and 25 April, 
2010 and on 7 May, 2010. Between 23 and 75 stan- 
dard stars (Landolt and SDSS) were used to establish the 
transformations to the standard photometric systems. 

2.4.2. Allegheny Lightcurve 

We obtained 44 nights of observations of GSC 03546- 
01452 with the 16-inch Keeler RCX-400 Meade telescope 
at the Allegheny Observatory, University of Pittsburgh. 
The observations span 18 months from May 2010 through 
November 2011 and were all made through a Cousins 
R filter onto an SBIG KAF-6303E/LE 2048x3072 pixel 
CCD with a pixel scale of 0.57" pixel -1 for a total FoV of 
19.5' by 29.2'. Exposure times ranged from 30-150 sec- 
onds depending on conditions, with a median exposure 
time of 75 seconds. 

Standard bias and dark subtraction and Hatfield cal- 
ibrations were preformed on the raw data, and the im- 
ages were astrom etrically calibrated to t he 2MASS Point- 
Source Catalog (|Skrutskie et al.ll2006D . Typical seeing 
conditions were between 3"-4" with a median seeing of 
3.5". Photometry was accomplished using a 10-pixcl 
(5.7") circular aperture and subtractring the estimated 
sky flux based on a 15-20 pixel (8.6"-11.4") radius sky 
annulus. Typical 1-er uncertainties were 5 millimag. Rel- 
ative photometry was calculated relative to the aver- 
age counts from two reference stars in the image with 
J2000 coordinates: (1): 19:11:34.733 +48:34:54.77; and 
(2) 19:11:56.335 +48:22:55.38. 



2.4.3. SuperWASP Lightcurve 

The SuperWASP survey (|Pollacco et al.l 120061 ) is a 
wide-angle transiting planet survey that monitors the 
brightness of millions of stars. The survey uses a vi- 
sual broad band filter that covers the wavelength range 
from 400-700nm. For details on reduction te c hniqu es and 
survey design please consult iPollacco et al.l (J2006D . The 
SuperWASP photometry for GSC 03546-01452 consists 
of 8823 points spanning just over four years from May 
2004 to Aug 2008. 

2.5. High Spatial Resolution Imaging 

Two high spatial resolution imaging campaigns of GSC 
03546-01452 were undertaken in order to help rule out 
false positives and to look for hierarchical structure. In 
particular, these observations were conducted to search 
for visual companions at large separations that could in- 
fluence our spectroscopic results, or that might be bound 
tertiary companions to GSC 03546-01452. The adaptive 
optics imaging and lucky imaging runs are complimen- 
tary with the adaptive optics focusing near the star, and 
lucky imaging covering out to larger distances from the 
star. 



? 3 

< 

" 4 
<0 

£ 5 

o 



GSC 03546 

NIRC2 AO, K' -filter 

August 31, 2011 



0.5 1.0 1.5 2.0 

angular separation (arcseconds) 

Figure 2. Contrast generated by NIRC2 AO observations of GSC 
03546-01452 . Off-axis sources with relative brightness AK' = 7 
are ruled out at 10-ct for angular separations beyond at 0.5". 

2.5.1. Adaptive Optics Imaging 

We also acquired adaptive optics (AO) observations 
of GSC 03546-01452 to assess its multiplicity at wide 
separations. Images were obtained on 31 Aug., 2011 
UT using NIRC2 (PI: Keith Matthews ) at the 10m 
Keck II telescope (jWizinowich et al.|[2000T) . GSC 03546- 
01452 (V=11.7) is sufficiently bright to serve as its own 
natural guide star. Our observations consist of 9 dithered 
images (10 coadds per frame, 0.5s per coadd) taken with 
the K' filter (A c = 2.12 fim). We used NIRC2's narrow 
camera setting, which has a plate scale of 10 mas pixel -1 , 
to provide fine spatial sampling of the instrument point- 
spread function. Raw frames were processed by clean- 
ing hot pixels, flat-fielding, subtracting background noise 
from the sky and instrument optics, and aligning and 
coadding the results. No off-axis sources were noticed in 
individual frames or the final processed image. Figure [2] 
shows the contrast levels generated by the observations. 
Our diffraction-limited images rule out the presence of 
companions down to AK' = 5.7, 7.1, 7.5 mags at angular 
separations of 0.25", 0.5", and 1.0" respectively. 

2.5.2. Lucky Imaging 

Lucky imaging (LI) involves acquiring large numbers 
of observations with a rapid cadence. A subset of these 
images are then shifted and stacked in order to produce 
nearly-diffraction-limited images. GSC 03546-01452 was 
obser yed over 3 nights se parated by 2 years using Fast- 
Cam (|Oscoz et al. Il2008( ) on the 1.5 m TCS telescope at 
Observatorio del Teide in Spain. The LI frames were ac- 
quired on 9 Oct., 2010; 25 Aug., 2011 and 11 Sep., 2012 
in the / band and spanning ~ 21 x 21 arcsec 2 on sky. 

The 11 Oct., 2010 observing run had 100,000 frames 
with 50 msecs per frame, the 25 Aug., 2011 run had 
60,000 frames with 60 msecs per frame, and the 11 
Sept., 2012 run had 75,000 frames at 60 msecs per frame. 
The data were processed using a custom IDL software 
pipeline. After identifying corrupted frames due to cos- 
mic rays, electronic glitches, etc., the remaining frames 
are bias corrected and flat fielded. 

Lucky image selection is applied using a variety of 
selection thresholds based on the brightest pixel (BP) 
method. This method involves selecting frames for the 
final combined LI image based on the BP in that frame. 
The selected BP must be below a specified brightness 
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Figure 3. GSC 03546-01452 was observed on 9 Oct., 2010, 25 
Aug., 2011, and 11 Sep., 2012 with LI from the Fastcam imager. 
The mosaic shows the final LI combined images from each run 
using different thresholds. LI detected a previously known visual 
companion 7.7" away from GSC 03546-01452 (marked with an ar- 
row) in the 9 Oct., 2010 and 11 Sep., 2012 runs. The companion 
was not detected on the 25 Aug., 2011 night because it was not in 
the field of view of the camera. This visual companion has AI m 
7.9 magnitudes. Further analysis showed that the co mpan ion is in 
fact a background star. For more details, see Section 1.3.51 

threshold to avoid selecting cosmic rays or other non- 
speckle features. As a further check, the BP must be 
consistent with the expected energy distribution from a 
diffraction speckle under the assumption of a diffraction- 
limited PSF. The BP's of each frame are then sorted from 
brightest to faintest, and the best X% are then shifted 
and added to generate a final combined image. Several 
values of X, the LI threshold, are tried until the best 
combined image is generated. The combined images us- 
ing different LI thresholds are shown for each of the three 
observing runs in Figure |3l The 11 Sept., 2012 run has 
a better quality image than the other two runs because 
of an instrument upgraded to a higher sensitivity CCD. 
The best LI threshold was determined to be 50%, which 
is what we used for the rest of the analysis. This results 
in a total exposure time for each of the final combined 
images from the three runs in order of date 2500, 1800, 
and 2250 s respectively. 

In Figure 0] we show the 3-c detectability curves com- 
puted out to 8" from GSC 03546-01452. W e follow the 
same procedure as in iFemenia et al. I (|2011l) to compute 
these curves: at a given angular distance p from GSC 
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Figure 4. Top: 3-<r detectability curves for LI from the FastCam 
imager. There was significant variation in the curves over the three 
runs. The best run can detect sources down to AI = 8.5 beyond 
2". No source was detected within 7" of GSC 03546-01452. Bot- 
tom: Conversion of 3-<r detectability curves for the three observing 
nights into mass sensitivities using empi rical Mass-Luminosity re- 
lationships in the literature. See Section 13.51 



03546-01452 we identify all possible sets of small boxes 
of a size larger but comparable to the FWHM of the PSF 
(i.e. 5 x 5 pixel boxes) . Only regions of the image showing 
structures easily recognizable as spikes due to diffraction 
of the telescope spider and/or artifacts on the read-out 
of the detector are ignored. For each of the valid boxes 
on the arc at angular distance p the standard deviation 
of the image pixels within the 5-pixel boxes is computed. 
The value assigned to the 3-ct detectability curve at p is 
3 times the mean value from the standard deviations of 
all the eligible boxes at p. This procedure, using each 
of the LI % thresholding values provides a detectability 
curve, while the envelope of all the family of curves for 
a given night yields the best possible detectability curve 
to be extracted from the whole data set. 

Although there are no companions detected within 7", 
we do find a possible companion slightly outside this re- 
gion at a separation of 7.7". This companion is signifi- 
cantly dimmer than GSC 03546-01452 with a AI « 7.9 
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magnitudes. This object was only detected in the 11 
Oct., 2010 and 11 Sept., 2012 campaigns. The visual 
companion was not detected in the 25 Aug., 2012 cam- 
paign due to the orientation of the camera on the sky. 
The possible companion can be seen in Figure [3] A dis- 
cussion of the likelihood that this object is actually a 
bound companion to GSC 03546-01452 can be found in 
Section [3.51 

3. RESULTS 
3.1. Host Star Characterization 



3.1.1. Spectroscopic An 

The ARCES moderate resolution spectrum (Section 
I2.3|) was analyzed by two independent analysis pipelines. 
We refer to these pipeline results as the "IAC" (Insti- 
tute de Astroffsica de Canarias) and "BPG" (Brazilian 
Participation Group) results. Briefly, both methods are 
based on the principals of Fe I and Fe II excitation and 
ionization equilibria. Both techniques e mploy the 2002 
version of the MOOG code (jSnedenlll973l) . but use differ- 
ent line lists, model atmosphere grids, equivalent width 
measurements, and convergence criteria. For a com- 
plete description of the analysis process, please refer to 
IWisniewski et all (|2012f ). 

Our IAC analysis used 204 Fe I lines and 20 Fe II 
lines, and resulted in the following parameters: T e g = 
5502±100, \ogg = 4.21±0.58, and [Fe/H] = +0.31±0.16. 
The BPG analysis used 61 Fe I lines and 6 Fe II lines re- 
sulting in these parameters: T c ff = 5652 ± 75, logg = 
4.46 ± 0.16, and [Fc/H] = +0.44 ± 0.10. Since the val- 
ues for both methods are consistent to within 1-ct they 
were combined using a weighted average as described in 
IWisniewski et al. (2012). This calculation resulted in the 
final spectroscopic parameter values of T c ff = 5598 ± 63, 
log 5 = 4.44 ± 0.17, and [Fe/H] = +0.40 ± 0.09. These 
values are recorded in Table and are used for the re- 
maining analysis. 

From these value s of T c g, logg , and [F e/H] and their er- 
rors, and using the lTorres et al.l (|2010l) relations, includ- 
ing intrinsic scatter and the uncertainties in the polyno- 
mial coefficients and covarianccs, we find that the mean 
and RMS of the mass and radius of the host arc M* — 
1.11 + 0.11 M Q and i?» = 1.06 + 0.23 R Q The median and 
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68 percent confidence intervals are M, 

and i?* = 1.03to:i9 R@- 

These spectroscopic values also allows us to transform 
the detection curves from the left panel of Figure 0] into 
upper limits on the mass of a possible stellar compan- 
ion undetected at the 3-ct level. We start by taking 
the spect roscopic T^a from Tabic [3l and using the ta- 
bles from[Mamaick] (J2010I ). we derive absolute V and I- 
band magnitudes. We can apply this technique because 
GSC 03546-01452 has been identified as a main-sequence 
star by the logg measurement. The Mj and the 3-tr dc- 
tectability curves allows the construction of the M/ vs. 
angular distance, p from the central star. This curve 
provides an absolute upper I-band limit to any com- 
panion, which would also have to be a main-sequence 
star. From the My vs. p curves we ca n use empiri 
cal Mass Luminosity Relat ionships (MLR) ()Henrv et al 
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Figure 5. The spectral energy diagram (SED) for GSC 03546- 
01452. The black line shows the best fit NextGen model compared 
to the observed flux in the photometric passbands from Tablc|3](rcd 
crosses). The horizontal bars are the approximate passband width 
for each filter, while the vertical bar is the error in the flux. The 
blue circles are the expected flux values from the model. The best 
fit model is fixed to the spcctroscopically determined T B ff, logg, 
and [Fc/H] and allowed Ay to float. The optical and infrared fit 
well, but there is a significant excess of flux in the GALEX NUV 
passband. 

the right panel of Figure |4j 

3.1.2. Photometric Analysis 

We corroborated the spectroscopic results with a series 
of photometri c tests, including determining the RPM-J 
statistic from iCollier Cameron et al.l (|2007l) . For GSC 
03546-01452, the RPM-J value was 2.31 and its (J-H) 
2MASS color was 0.387, results consistent with GSG 
03546-01452 being a dwarf. 

We also constructed a spectral energy diagram (SED) 
of GSC 03546-01452 using the GALEX NUV filter 
(|Morrissev et al.l 12007ft ; Johnson B and SP SS u' ,q' ,r' ,i' 
obser vations from HAO; 2MASS J, H, and K dCutri et al 
20f)l: and WISE1, WISE2, WISE3 (jCutri fc et al 



Xia fe Full2010D to derive the upper mass limit. These 
constraints arc plotted for all three observation nights in 



2012aj ). A summary of the photometric values are pre- 
sented in Ta ble EH We used the Ne xtGen model atmo- 
sphere grid (jHauschildt et all 11999!) to construct theo- 
retical SEDs. These models were fixed by the spec- 
troscopic values for T^Eff, log.g, and [Fe/H] described 
in Section 13.1.11 the reddening was constrained by the 
iSchlegel et al.l (|1998| ) dust maps in the direction of (l,b) 
= (79?339029, 16?853379) to be less than Ay = 0.213. 
The resulting fit to the photometry is shown in Figure 
[5j The best fitting model has a reduced \ 2 of 3.91, a 
negligible reddening of Ay = 0.1 ± 0.1, and a distance 
of 219 ± 21 pc. There appears to be a slight UV excess, 
which may be the result of modest stellar activity. 

3.1.3. Evolutionary Status and and Galactic Population 

Using the spectroscopic host star parameters and the 
mass and radius derived from the iTorres et alJ (|2010l ) 
relations, GSC 03546-01452 can be placed onto an 
evolutionary track. We use the Yonsei- Yale ( "Y2" ) 
model tracks from iDemarque et al.l (|2004l ) (and refer- 
ences therein), and select the track corresponding to 1.11 
M Q with a [Fe/H] = +0.40. Figure [6] shows this track 
with stellar ages marked as blue dots on the track. The 
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Table 3 

Host Star Properties GSC 03546-01452 



Parameter 



Value 



\-a Uncertainty 



Source 



GSC 1.1 Name 
GSC2.3 Name 
KIC Name 
2MASS Name 
a (J2000) [deg] 
8 (J2000) [deg] 
fi a [mas yr _1 ] 
fig [mas yr _1 ] 

T eff [K] 

log g [cgs] 
[Fe/H] [dex] 
Vmic [km s _1 ] 
Vrot sin i [km s" 

^systemic [Km S 

A v [mag] 
Distance [pc] 
M, [M Q ] 
R* [R©] 



3546-01452 

N2EA000110 

11022130 

J19113252+4830436 

287.885735 

48.512117 

-28.6 

-8.7 

5598 

4.44 

+0.40 

0.38 

<9.0 

-13.799 

0.1 

219 

1.11 

1.06 



0.26" 

0.25" 

2.3 

1.7 

63 

0.17 

0.09 

0.17 

0.034 

0.1 

21 

0.11 

0.23 



Laskeretal 

^askej^et^al 

13rowri ii et_al 

Cutrietal. 

Laskeretal 

Laskeretal 

Zachariasetal 

Zachariasetal 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 

This work 




galNUV 

B 

V 

Ro 

lc 

SDSS u' 

SDSS g' 

SDSS r' 

SDSS V 

2MASS J 

2MASS H 

2MASS K s 

WISE1 

WISE2 

WISE3 



18.320 
12.593 
11.799 
11.335 
10.967 
13.483 
12.156 
11.560 
11.400 
10.46 
10.07 
10.010 
12.650 
13.344 
15.181 



0.100 
0.020 
0.010 
0.012 
0.012 
0.020 
0.012 
0.010 
0.010 
0.030 
0.030 
0.03 
0.023 
0.020 
0.032 



Lorri ssey i 
HAO (This woTKJ 
Derived from HAO 
Derived from HAO 
Derived from HAO 
HAO (This work) a 
HAO (This work) a 
HAO (This work) a 
HAO (This work) a 
Cutri et al. (2003*1 
Cutri et al. (2003 
Cutri et al. (2003 ) 

Cutri&etaL 

^jutri_j^^t_aL 

Cutri fc et al. 



a These passbands were used to generate SED in Figure [5] 

dashed lines show tracks for the same metallicity, but for 
±O.11M , which is the 1-tr uncertainty on our mass es- 
timate. The gray area shows the region that this family 
of tracks occupies. The red point shows GSC 03546- 
01452 with 1-er error bars. From the error bars in logg 
and T c ff in Figure [6] we can constrain the age of GSC 
03546-01452 to being less than approximately 6 Gyr old. 
We can determine the Galactic population member- 
ship of GSC 03546-01452 by using the absolute systemic 
RV=-13.799 ± 0.128 km s" 1 (where the 0.128 km s" 1 
is a conservative estimate of the systematic error from 
the comparison wit h RV standard stars) , the proper mo- 
tions from UCAC4 (jZacharias et alJl2012l ) (fi = -28.6 ± 
2.3, — 8.7 ± 1.7 mas yr _1 ), and the distance from Tablc|3] 
(219 ± 21 pc). We find (U, V, W) = (24.0 ± 2.5, -10.1 ± 
1.3, 25.6 ± 3.1) km s _1 . This velocity is consistent with 
m embership in t he thi n disk according to the criteria 
of iBensbv et al.l (|2003f ). The space motion velocities 
were determined using a modification of the IDL rou- 
tine GALJJVW, which is ultim ately based on the method 
Uohnson fc Soderbloml ([1987D . We adopt the correction 
for the Sun's mot ion with respect t o the Local Stan- 
dard of Rest from iCoskunoglu et al.l (|2011[ ) , and choose 
a right-handed coordinate system such that positive U is 
toward the Galactic Center. 

3.2. Companion Orbital Analysis 
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Figure 6. The black line represents the best fit Yonsci-Yale 
evolutionary track for a 1.1 Mq star with a [Fe/H] of +0.40 
(Dcmarque et al. 2004). The gray area denotes the \-a deviation 
from that track. Ages along the track are denoted by blue dots 
for 1.0, 4.7, 6.4, 8.0, and 9.0 Gyr. The red cross denotes the loca- 
tion of GSC 03546-01452 on this diagram with the spectroscopic 
uncertainties. 

The orbital parameters of MARVELS-6b were derived 
using the radial velocity data from both the MARVELS 
and the SARG spectrogr ap hs. T he orbital fit was cre - 
ated using the EXOFASTE3 code (jEastman et al.ll2013l) . 

27 http://astroutils.astronomy.ohio-state.edu/exofast/ 
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Table 4 

Properties for MARVELS-6b 



Parameter 



Units 



Value 



T c BJDtdb - 2450000 

P Period (days) 

e eccentricity 

u) argument of periastron (radians) 

K RV semi-amplitude (m/s) 

7TJVG TNG Systemic velocity (m/s) . . . 

dv/dt RV slope (m/s/day) 

"YAPO APO Systemic velocity (m/s) . . . 

e cos(cj) 

e sin(u) 

T P BJD TDB - 2450000 

Mi, sin i Minimum Mass (Mj up ) 



5023.377; 



-o.io 

-0.099 



47 C02Q+ ' 0063 
tl •°»^»_0.0062 

1442+ - 0078 
u- 1 ** -0.0073 

,+0.064 

'-0.061 



1. 



1611 



-12 
-13 



724 ± 15 



0.180 ±0.053 

1058 
-0.0597 
0.1312 
5025.81 
31.7 



+ 17 
16 

0.0071 
0.0070 
+0.0093 
0.0092 
+0.46 
-0.44 
2.0 
12.(1 



which allowed us to determination of the usual Keplerian 
orbital elements shown in Table @] 

The first step in the process was to perform indepen- 
dent fits to both the MARVELS and SARG data points. 
Both data sets had error bars that were too large based 
on the x 2 of their fits, so we re-scaled them to force the 
probability that the \ 2 is greater than or equal to the 
measured x 2 , i.e., P(x 2 ), = 0.5. The MARVELS error 
bars were scaled by a factor of 0.579, and the SARG er- 
ror bars were scaled by a factor of 0.222. The error bars 
in Table [T] include these scaling factors. 

Once the error bars were scaled, the two data sets were 
combined. The final combined radial velocity curve was 
then run through EXOFAST which is a Markov Chain 
Monte Carlo (MCMC) based IDL software suite. The 
basic algorithms of this software follow the work of lFordl 
(|2006| ) . The results of the combined fit arc listed in Table 
21 and the best-fit model is shown in Figures [7] and [8] 

The radial velocity measurements from both MAR- 
VELS and SARG are relative radial velocities. As part 
of the fitting process, EXOFAST fits and subtracts arbi- 
trary zero points for each data set simultaneously. These 
zero points are recorded in Table |U These offsets are 
based on arbitrary zero-points, and should not be con- 
fused with the true systemic velocity of the host system 
-13.799 ±0.128 km s" 1 . 

The systemic velocity of GSC 03546-01452 was derived 
by taking the orbital solution from EXOFAST for the 
combined relative radial velocity data set and applying 
it to the SARG/TNG absolute velocity data set (Table 
[21). This calculation was performed by using the MPFIT 
Levenberg-Marq uardt non-linear le ast squares fitting al- 
gorithm in IDL (|Markwardtll2009[ ). A Keplerian model 
was fit to the SARG/TNG radial velocity points, holding 
P, e, u>, K, and Tp fixed while allowing 7 to float; the 
results are presented in Table El 

3.3. Companion Mass 

3.3.1. Mass Functions of Secondary 

The mass function (A4b) is the only property of a 
single-lined radial velocity companion that we can de- 
rive that is independent of the properties of the primary. 
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Figure 7. All the radial velocity points from both the MARVELS 
Survey (blue squares) and TNG/SARG follow-up spectroscopy (red 
circles). The final best fit for a M^sini = 31.7 ± 2.0 Mj up com- 
panion to GSC 03546-01452 (Table [4j is shown with the associ- 
ated O-C diagram. The observation data points start at BJD = 
2454956.916 (5 May, 2009). There is a slight residual linear slope 
(0.180 m s _1 day -1 ) after removing the companion orbit. 

Using the MCMC chain from the joint RV fit, we can 
derive Mj, of the companion, 



M 



(AfhSini) 3 
= {M*+M b y 

= A" 3 (l-e 2 ) 3 / 2 



P 



2nG 
: (2.136 ±0.049) x KT 5 M , 



(1) 



where the uncertainty is essentially dominated by the 
uncertainty in K, such that OM b jM.\ i ~ 3(ctk/K) — 
3 x 0.8% ~ 2.4%. 

3.3.2. Minimum mass and mass ratio 

In order to determine the mass or mass ratio of the 
secondary, we must estimate the mass of the primary, as 
well as the inclination of the secondary. 

To estimate the mass and radius of the primary, we 
proceed as follows. For each link in the MCMC chain 
from the joint RV fit, we draw a value of T e g, log g, and 
[Fc/H] for the primary from Gaussian distributions, with 
me ans and dispe r sions given in Table [31 We then use 
the iTorres et al.l ([2010D relations to estimate the mass 
M* and radius i?* of the primary, including the intrinsic 
scatter in these relations. 

The minimum mass (i.e., M b if sin i = 1) and minimum 
mass ratio of the secondary are: 



M fe , min = 31.7 ± 2.0 M Jup 

= 0.0303 ± 0.0020 M 
q m in = 0.02733 ±0.00092 



(2) 



The uncertainties in these estimates are almost entirely 
explained by the uncertainties in the mass of the primary: 
o- Mb /M b ~ (2/3)(ov./M») = (2/3) x 9.7% ~ 6.5%, close 
to the uncertainty in M b mul above (6.4%), and cr q /q ~ 
(1/3)(Vm,/M,) - (1/3) 'x 9.7% - 3.2%, close to the 
uncertainty in q (3.4%). 

3.3.3. A posteriori distributions of the true mass 

The a posteriori distribution of the true mass of the 
companions given our measurements depends on our 
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Figure 8. The phased final RV curve including both the MAR- 
VELS Survey (blue squares) and the TNG/SARG observations 
(red circles). The final best fit for a M^sini = 31.7 ± 2.0 Mj up 
companion to GSC 03546-01452 (Table|4j, is shown with the asso- 
ciated O-C diagram. 

prior distribution for the mass of the companion, which 
is roughly equivalent our prior on the mass ratio. 

If we assume a prior that is uniform in the logarithm 
of the true mass of the companion, then the posterior 
distribution of cos i will be uniform. More generally, for 
other priors, cosi is not uniformly distributed. We adopt 
priors of the form: 

dN * /»> 

^ (Xq (3) 

where q is the mass ratio between the companion and the 
primary, and a = — 1 for the uniform logarithmic prior 
discussed above. To include this prior, for each link the 
MCMC chain we draw a value of cos i from a uniform dis- 
tribution, but then weight the resulting values of the de- 
rived parameters for that link (i.e., the companion mass 
m) by q a+l . 

For a > 0, the a posteriori distribution does not cov- 
erage. However, we can rule out equal mass ratio com- 
panions by the lack of a second set of spectral lines. We 
therefore assume that the mass ratio cannot be greater 
than unity; i.e., we give zero weight to inclinations such 
that q > 1. In doing so, we implicitly assume that the 
companion is not a stellar remnant. 

Figures O shows the resulting cumulative a posteri- 
ori distributions of the true mass of the companion, for 
a = — 1 (uniform logarithmic prior on q) , a = (uniform 
linear prior), and a = 1. For a = — 1 and a = 0, the 
inferred median masses are ~ 37 Mj up and ~ 45 Mj upi 
respectively. For these priors, we conclude that the com- 
panion has a mass below the hydrogen burning limit and 
thus is a bona fide brown dwarf at 90% and 72% con- 
fidence, respectively For a = 1, the median mass is 
172 Mj up or ~ 0.16 M Q , firmly within the stellar regime. 
Indeed, for this prior, there is only a ~ 32% probability 
that the companion is a brown dwarf. However, this con- 
clusion is sensitive to the precise form for our constraint 
that q < 1. With a more careful analysis it may be 
possible to rule out a wider range of companion masses 
based on the lack of evidence for a second set of spectral 
lines. Furthermore, there is little evidence that the mass 
function of companions to solar type stars is rising as 
steeply as a = 1 for minimu m masses around 30 Mj up 
(|Grether fc Lineweaverl 120061 ) . 

We conclude that the companion to GSC 03546- 
01452 is most likely a brown dwarf. However, we cannot 
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Figure 9. The cumulative a posteriori probability of the true 
mass of the secondary, for three different priors on the companion 
mass ratio, dN/dq oc q a , with a = —1 (solid), a = (dotted) 
and a = 1 (dashed). Also indicated are the mass of the primary 
with uncertainty and the hydrogen burning limit both denoted by 
shading. For a = — 1 and a = 0, the inferred median masses 
arc ~ 37 Mj up and ~ 45 Mj up , respectively. For these priors, 
we conclude that the companion has a mass below the hydrogen 
burning limit and thus is a bona fide brown dwarf at 90% and 72% 
confidence, respectively. For a = 1, the median mass is 172 Mj up 
or ~ 0.16 Mq, firmly within the stellar regime. However, it is 
unlikely that a > in this mass regime. 

definitively exclude that it is, in fact, a low-mass stellar 
companion seen at low inclination. 

3.4. Time-Series Photometric Analysis 
3.4.1. Summary of Datasets 

The SuperWASP photometric dataset for GSC 03546- 
01452 consists of 8823 points spanning just over four 
years from HJD' = 3128 to 4690 (where HJD' = HJD - 
2450000). The full, detrended SuperWASP dataset has 
a relatively high weighted RMS of 1.8% and exhibits evi- 
dence for systematics. The distribution of residuals from 
the weighted mean is asymmetric and non-Gaussian, 
showing long tails containing a much larger number of 
> 3-er outliers than would be expected for a normally- 
distributed population. We therefore cleaned the Super- 
WASP data as follows. We first add a trial systematic 
fractional error in quadrature to the photon noise un- 
certainties o~ sys . We then compute the error- weighted 
mean flux, and determine and reject the largest error- 
normalized outlier from the mean flux. We then recom- 
pute the mean flux, and scale the uncertainties by a con- 
stant factor r to force x 2 /dof = 1. We iterate this pro- 
cedure until no more > 4-er outliers remain. We then 
repeat the entire procedure to determine the value of 
a sys that results in a distribution of normalized residu- 
als that has the smallest RMS from the Gaussian expec- 
tation. Although 4-er is a slightly larger deviation than 
we would expect based on the final number of points, 
we adopt this conservative threshold in order to avoid 
removing a potential transit signal. We adopt r = 0.60 
and a sys = 0.011. The cleaned light curve has 8635 data 
points with an RMS of 1.5% and x 2 /dof = 1 (by design). 

The Allcngheny photometric dataset consists of 3990 
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Figure 10. Cleaned broad band relative photometry of GSC 
03546-01452 from SuperWASP (blue) and Allegheny (green) as a 
function of date. 

points spanning roughly a year and a half from 
HJD'=5341 to 5884. The weighted RMS of the raw light 
curve is 0.63%. There is mild evidence for systematics 
errors in this dataset, and thus we repeat the identical 
procedure as with the SuperWASP data. We choose a 
more conservative 6-cr cut, given the better precision of 
the Allenghcny observations, to prevent us from remov- 
ing real transit-like variability. We adopt r = 1.04 and 
a sys = 0.0028, with a final RMS of 0.4% from 3983 data 
points. 

While continuous variability will not be removed by 
our procedure, variability in the form of a small number 
of highly-discrepant point will be masked. This possibil- 
ity is particularly relevant for transit signatures. Given 
the typical uncertainties of ~ 0.4% for that data set, 
we do not expect our cleaning procedure to remove out- 
liers that differ by less than ~ 2.4% from the mean, 
roughly corresponding to the typical depth for a tran- 
sit of a ~ 1.5 Rjup companion given the primary radius 
of- 1.03 R Q . 

Finally, we combine all the relative photometry af- 
ter normalizing each individual data set by its mean 
weighted flux. Figure [TU] shows the combined data set, 
which consists of 12618 data points spanning roughly 7.6 
years from HJD' = 3128 to 5884. The weighted RMS 
is 0.66%. The resulting light curve is constant to within 
the uncertainties over the entire time span. Because the 
individual data sets sample disjoint times and are not 
relative common set of reference stars, our procedure of 
normalizing each data set by its mean flux will remove 
variations on the longest timescales. 

3.4.2. Search for Variability 

We ran a Lomb-Scargle periodogram on the full pho- 
tometric dataset, testing periods between 1 — 10 3 days. 
The resulting periodogram shown in Figure QT] does not 
exhibit any strong peaks. The maximum fitted ampli- 
tude over this range is ~ 0.1%. The inset presents the 
periodogram of the combined data set for periods within 
— 10 days of the period of the companion. The individ- 
ual periodograms for the SuperWASP and Allcngheny 



Figure 11. Lomb-Scargle periodogram of the combined photo- 
metric data. While a large number of peaks are visible, we do not 
regard any of them as significant. In particular, while there is a 
local peak at the period of the secondary (vertical red dotted line) , 
this peak arises from the SuperWASP data alone, and is not con- 
firmed by the Allegheny data. This result is demonstrated in the 
inset, which shows the combined periodogram near the period of 
the companion (grey solid), the periodogram of the SuperWASP 
data alone (blue dotted), and Allegheny alone (blue long dashed). 



data are also shown. While the SuperWASP data dis- 
play a local peak at the period of the companion, this 
feature is not corroborated by the Allengheny data. Al- 
though this diffrence could in principle arise from real 
variability that is present in the SuperWASP data, but 
not in the Allegheny data (i.e., evolving spots), given 
the higher precision of the Allcngheny data, we believe 
it is more likely that the peak in the SuperWASP peri- 
odogram is due to low-level systematics in the form of 
residual correlations on a range of time scales. 

Figure [12] shows the combined light curve, folded at the 
median period and time of conjuction of the companion, 
(P = 47.8929 ± 0.0063 and T c = 5023.377), and binned 
0.05 in phase. The weighted RMS of the binned light 
curve is ~ 0.041%. Although the variations are larger 
than expected from a constant light curve based on the 
uncertainties (x 2 /dof — 2.6), we again suspect that these 
are due to systematic errors in the relative photometry. 
In particular, the folded, binned Allegheny light curve 
shows a somewhat lower RMS of — 0.016%, and is more 
consistent with a constant flux with a x 2 /dof = 0.4. 

We conclude that there is no strong evidence for vari- 
ability of GSC 03546-01452 on any time scale we probe. 
We can robustly constrain the amplitude of any persis- 
tent periodic variability to be less than < 0.1% over 
timescales of < 3 years, and we can constrain the ampli- 
tude of persistent photometric variability at the period 
of the companion to < 0.02%. 

3.4.3. Excluding Transits of the Secondary 

For a uniform distribution in cos i (corresponding to a 
prior that is uniform in the logarithm of the mass), the 
secondary transit probability is — 1.8%. The duration 
for a central transit is — i?*P/(7ra) — 6.6 hours, and 
the depth for Jupiter-radius body is 5 — (Rb/R*) 2 = 
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Figure 12. Relative photometry folded at the period of sec- 
ondary, and binned 0.05 in phase. Phase zero corresponds to the 
expected time of conjunction (and so of transits for the appropratc 
inclinations) . Black points are the combined data, blue are Supcr- 
WASP, and green are Allegheny. 

1.0%(Rb/Rj U p) 2 , where Rb is the radius of the compan- 
ion, and we have adopted the median value for the pri- 
mary radius of i?» = 1.03 Rq. Thus, the expected S/N of 
a transit in the combined photometric dataset assuming 
uniform phase coverage is 



S/N 



N x/i i'R» 



1/2 



£-13 

a 



R\, 



R 



J up 



(4) 



where N = 12618 is the number of data points, and 
a ~ 0.66% is the RMS of the combined light curve. Thus 
we would expect to be able to robustly detect or exclude 
transits from companions with radii as as small as ~ 
0.9 Rjup at 10-cr if the transit phase is well covered by 
the photometric data. However, given the long period of 
the companion, uniform phase coverage is not necessarily 
a good approximation. 

Therefore, in order to account for the actual sam- 
pling of the photometric data, we perform a quantitative 
search for transit signals using a modified Monte Carlo 
method. This meth od is similar to that described in 
iFleming et al.1 (J2012I ) we briefly review the details here. 
We start with the distributions of the relevant radial ve- 
locity fit parameters: the period P, scmiamplitude K, 
eccentricity e, and time of conjuction Tq- These are ob- 
tained from the MCMC chain determined by fitting of 
the joint radial velocity data set as described in 13.21 For 
each link in this MCMC chain, we also draw values for 
T ff, logg, and [Fc/H] for the primary from Gaussian 
distributions, with mea ns and dispe r sions given in Ta- 
ble [31 We then use the iTorres et al.1 (|2010l ) relations to 
estimate the mass M* and radius i?* of the primary, in- 
cluding the intrinsic scatter in these relations. We then 
draw a value of cos i from a uniform distribution, assum- 
ing a prior on the companion mass Mb that is uniform 
in log Mb. The values of P, K, e, M*, i?» and i are 
then used to determine the secondary mass Mb, semi- 
major axis a, and impact parameter of the secondary 
orbit b = a cos i/R*. From the impact parameter, we 



determine if the companion transits (b < 1), and, if so, 
we determi ne the properties of the light curve using the 
routines of Mandel & Agol (2002) for a given compan- 
ion radius Rb- We assume quadratic limb darkening, 
adopting coeffi cients appropriate for the R band from 
iClaret fc Hauschild t (2003), assuming solar metallicity 
and the central values of T e g and logp listed in Table 
[3l Finally, we fit the predicted transit light curve to the 
combined photometric data, and compute the A% 2 be- 
tween the constant flux fit and the transit model. We 
perform this calculation for every step in the Markov 
Chain, and finally repeat this procedure for a range of 
companion radii. 

We find a best-fit transit model has Ax 2 = —9.0 rela- 
tive to a constant fit. We do not consider this result to be 
significant, as we find a larger improvement A^ 2 = —28.1 
when we consider signals with the same ephemeris and 
shape as transits, but corresp onding to positive devia- 
tions (i.e., "anti-transits", see I Burke et al.1 (120060 ). Al- 
though formally statistically significant, both the transit 
and anti-transit signals are likely caused by residual sys- 
tematics in the cleaned photometric data set, and thus 
we conclude there is no evidence for a transit signal in 
the photometric data. 

This procedure can be used to determine the confi- 
dence with which we can rule out transits of a compan- 
ion with a given radius. Specifically, the confidence with 
which we can exclude a companion with a given Rb is sim- 
ply the fraction of the steps in the Markov Chain where 
the companion transits, which produces a transit with a 
Ax 2 relative to the constant fit that is greater than some 
value of Ax 2 - We consider three different thresholds of 
Ax 2 = 25, 49, and 100. The resulting cumulative prob- 
ability distributions are shown in Figure 1131 Given the 
values of Ax 2 found for anti-transits, we conservatively 
consider thresholds of A^ 2 > 49 to be robust. For this 
value, we can exclude companions with Rb > 0.7 Rj up 
with 50% confidence (i.e., for 50% of the trials), and com- 
panions with Rb > 1.2 Rj with 90% confidence. How- 
ever, there is a long tail toward large companion radii, 
arising from the conflation of the long period of the com- 
panion, uncertainties in the emphemeris, and imperfect 
phase coverage. Therefore, we are unable to exclude 
transits at > 95% confidence even for Rb > 1.5. For the 
minim um companion mass of ~ 32 Mj up , iBaraffe et al.l 
(|2003| ) predict radii of ~ 0.9 — 1 Rj up for ages between 
0.5 and 5 Gyr. Therefore, we cannot definitively exclude 
the possibility that the companion transits. 

As discussed above, our 6-<r clipping procedure would 
in principle remove > 2.4% transit signatures from the 
Allegheny dataset arising from a large Rb > 1.5 Rj U p 
companion. We therefore repeated the transit search on 
the Allegheny dataset with 10-cr clipping, but again find 
no significant transits. 

3.5. LI Tertiary Companion 

As discussed in Section 12.5.21 a possible tertiary com- 
panion was detected at 7.7" from GSC 03546-01452 (see 
Figured]). The possible tertiary companion has AI « 7.9 
magnitudes, and at the measured distance to the host 
star in Table [3] places the companion at a projected dis- 
tance of 1686 AU from the host star. 

This possi ble tertiary comp anion can be found in both 
the 2MASS ()Cutri et al.ll2003l ) and Kepler Input Catalog 
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Figure 13. Probability that transits of a companion are excluded 
at levels of Ax = 25, 49, 100 based on the analysis of the combined 
SupcrWASP and Allegheny photometric data sets, as a function 
of the radius of the companion. The gray region represents the 
95 percent confidence level. Only the Ax 2 = 25 model crosses 
into the gray region, thus we cannot rule out a possible transiting 
companion even with an unphysically large radii R\, > 1.5 Rj U p at 
the 95% level. 



(KIC; iBrownet al]|201ll) . The KIC ids of the primary 
and possible tertiary are 11022130 and 11022139 respec- 
tively. The first blow to the hypothesis that this is a 
bound system is found in the KIC itself. The KIC gives 
the following stellar parameter estimates for the possi- 
ble tertiary T cff = 5919, log 5 = 4.361, and A v = 0.515. 
Comparing these values to the spectroscopic values de- 
termined for GSC 03546-01452 in Table [3] it is clear that 
these two objects cannot be bound, and have a AI of « 
7.9 magnitudes; because two G dwarfs should have ap- 
proximately the same magnitude. However, the tertiary 
is faint (J = 16.272), and the KIC spectroscopic determi- 
nations are photometric in nature, so we considered the 
possibility that the KIC parameters were incorrect. 

GSC 03546-01452 has a measured prope r motion from 
the UCAC4 catalog (jZacharias et al.ll2012t) of \x a = -28.6 
± 2.3 mas yr _1 and fis = -8-7 ± 1.7 mas yr _1 (see Ta- 
ble [3J . This motion is significant given the plate scale of 
the LI images ~ 42.56 mas pixel -1 , which, given our 2 
year baseline for observations, should be sufficient to test 
whether this is a bound system. Table [3] gives the differ- 
ence in position between the GSC 03546-01452 and the 
possible tertiary in RAcos(DEC) and DEC for both LI 
observations. The two points are distinct to the 2-cr level. 
We also have observations of both sources in 2MASS, 
although the astrometric position measurement is sig- 
nificantly worse, the 12 year baseline between 2MASS 
and the 2012 observation compensates for the lower ac- 
curacy. These points are also distinct to the 2-tr level. 
Given both the astrometric and photometric data, the 
hypothesis that this is a bound system is very unlikely. 

4. DISCUSSION 

Based on the derived spectroscopic parameters for the 
host star and the results discussed in Section [3.3.31 it is 
highly likely that MARVELS-6b is a BD companion to 
GSC 03546-01452 assuming a logarithmic or linear prior 
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Difference in Position of GSC 03546-01452 versus Possible Tertiary 

Year Observed Aocos(<5) Error A<5 Error Source 
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on the distribution of mass ratios. Working with the 
minimum mass for MARVELS-6b of 31.7 ± 2.0 M Jup , we 
can comp are it to current distribution of masses in the 
literature. iGrether fe Lineweaverl (J2006D investigated the 
mass function of both stellar and planetary companions 
and found that they both decreased by a substantial frac- 
tion as they approached the minimum of the BD desert. 
In the case of the stellar mass function, it dropped by 
two orders of magnitude from the 1 Mq to the middle of 
the BD desert. On the planetary mass side, it rises one 
order of magnitude as one progress away from the BD 
desert and towards lower mass planets. 

MARVELS-6b is interesting because its minimum 
mass, 31.7 ± 2 .0 Mj nr , lies at the very botto m of the two 
distributions. IGrether fe Lineweaverl (|2006[ ) define their 
mass functions as number of companions per unit interval 
in log mass, and find the minimum of the distributions 
to be at 31~^^Mj up . In order to better understand this 
minimum, it is necessary to increase the number of ob- 
jects (or place tighter constraints by non-detections) in 
this region of the mass function, and MARVELS-6b takes 
us one step closer to that goal. 

MARVELS-6b stands out in respect to other brown 
dwarfs in two significant ways. First, it has a relatively 
low eccentricity, 0.1442, given its moderately long pe- 
riod of si 47 days. To place this point in perspective, 
MARVELS-6b is plotted against a catalog of BDs from 
iMa fc Gel (|2013bf ) in Figure [H Two aspects of this plot 
are particularly interesting: First, MARVELS-6b is in 
an undcrdense region in period space. In the inset of 
Figure [14] there is a clear rise in the envelope of eccen- 
tricity as one moves to larger periods until one reaches 
roughly 200 days. MARVELS-6b is either on or below 
the bottom edge of this envelope. There is the possi- 
bility that this effect is due to observational bias; for 
low signal to noise ratios (SNR), one often measures a 
non-zero eccentricity even when the true eccentricity is 
zero. For a given radial velocity precision and compan- 
ion mass, longer periods will have lower SNR (since the 
semi-amplitude is smaller) , and so the lack of low eccen- 
tricity companions at long periods could be due to this 
effect. In any case, there are so few BD in this plot, that 
it is not clear whether MARVELS-6b's low eccentricity is 
significant, so further comment on this will have to wait 
until more data is available. 

This BD companion is also an outlier in metallicity 
space. In Figu re fi"5l MARyELS-6 b is again plotted with 
the BD from iMafc Gel (|2013bf ). GSC 03546-01452's 
metallicity is significantly super solar, and MARVELS- 
6b exists in an empty region of eccentricity-mctallicity 
space. As in the previous figure, the plot is merely sug- 
gestive of interesting astrophysics, but there are not a 
sufficient number of BDs to know whether MARVELS- 
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Figure 14. Orbital period and eccentricity of MARVELS-6b com- 
pared to a subset of literature brown dwarfs found through tran- 
sit or radial velocity that occupy the period range of to 100 
days. T he literature brown dwarfs are from the catalog of lMa fc Gel 
12013bf l. MARVELS-6b occupies an empty part of orbital period 
/ eccentricty plane. Inset: The full brown dwarf catalog extending 
out to much larger periods. 
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Figure 15. The host metallicity and orbital eccentricity of 
MARVEL S-6b comp a red to the literature brown dwarfs from the 
catalog of IMa fc Gel <2013bn MARVELS-6b has a relatively low 
eccentricity and its host star has a quite high [Fe/H] . 

6b is a significant deviation. MARVELS-6b is helping to 
fill in the BD parameter space, and with more data its 
location may become more astrophysically significant. 

5. CONCLUSION 

In this paper, we report the discovery of a BD compan- 
ion to GSC 03546-01452 with a minimum mass of 31.7 
± 2.0 M Jup . This BD, designated MARVELS-6b, has a 
moderately long period of 47.8929JIo 0062 days with a low 
eccentricty of 0.1442to;oo73- We have analyzed moderate 
resolution spectroscopy of the host star and have deter- 
mined the following properties: T e g = 5598 ± 63, log g = 
4.44 ±0.17, and [Fe/H] = +0.40 ±0.09. From these mea- 
surements we find that GSC 03546-01452 has a mass and 
radius of Ah = 1.11±0.11M Q and i?* = 1.06±0.23 R Q . 
This result combined with photometry, indicates the host 
star is a G dwarf at 219 ± 21 pc from the Sun with an age 
less than approximately 6 Gyr based on an evolutionary 
track analysis. 

Due to its moderately long period, MARVELS-6b has 
a transit probability for a uniform distribution of cosi 
of only 1.8%. Although we have roughly 13,000 pho- 
tometric data points, we cannot conclusively rule out a 
transit. In the Keck AO imaging no visual companions 
were found. However, in the LI a previously known com- 



panion at 7.7" from the host star was detected. This 
visual companion appears in both the 2MASS and KIC 
catalogs, and was shown not be a physical companion 
based upon photometry and astrometry, which is unex- 
pected given that many of the previous BD found by the 
MARVELS survey did have tertiary companions. 

Finally, we found that minimum mass of MARVEL S- 
6b exists at the minimum of the mass functions of close 
(orbital period < 5 yr) stellar and planetary compan- 
ions to stars, making this a rare object even compared 
to other BDs. It also exists in an undcrdense region 
in both period/eccentricity and metallicity/eccentricity 
space. This ultimately furthers the goal of this series of 
papers, to help fill out the low-mass companion phase 
space, which will ultimately help us understand these 
intriguing objects. 

The SuperWASP and Allegheny lightcurve data for 
GSC 03546-01452 , along with the APO spectroscopic 
data will be made available through the Vizier/ CD SP 8 ! 
catalog service. The SuperWASP and Allegheny 
lightcurve data will also be available with the online edi- 
tion of this article. 
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